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A parameter sensitivity study was conducted at the NASA Glenn Research Center's Icing Research Tunnel (IRT) 
using a 36 inch chord (0.91 m) NACA-0012 airfoil. The objective of this preliminary work was to investigate the 
feasibility of using ice shape feature changes to define requirements for the simulation and measurement of SLD 
icing conditions. It was desired to identify the minimum change (threshold) in a parameter value, which yielded an 
observable change in the ice shape. Liquid Water Content (LWC), drop size distribution (MVD), and tunnel static 
temperature were varied about a nominal value, and the effects of these parameter changes on the resulting ice 
shapes were documented. The resulting differences in ice shapes were compared on the basis of qualitative and 
quantitative criteria (e.g. - mass, ice horn thickness, ice horn angle, icing limits, and iced area). This paper will 
provide a description of the experimental method, present selected experimental results, and conclude with an 
evaluation of these results, followed by a discussion of recommendations for future research. 


Nomenclature 


AO A 

angle of attack 

IRT 

Icing Research Tunnel 

UVC 

liquid water content 

MVD 

median volume diameter 

SLD 

super-cooled large droplets 

T, 

total temperature 

T s 

static temperature 

V 

airspeed 

WBS 

work breakdown structure 


I. Introduction 

How good does it have to be? That is a question which has been asked repeatedly of scientists and engineers 
developing new capabilities for simulation and measurement of SLD icing conditions. The recent emphasis by 
NASA and other international partners on extending the current capabilities of engineering tools (icing tunnels and 
tankers, ice accretion computer codes) to include Super-cooled Large Droplets (SLD) has once again focused 
renewed attention on this question. 

NASA and it's international partners developed an SLD Technology Roadmap to guide the development of SLD 
Engineering Tools to support certification in SLD icing conditions. This SLD Technology Roadmap identified the 
key technology elements and supporting research investments needed to develop SLD engineering tools. 1 The SLD 
Technology Roadmap was later translated into a project plan having a schedule and Work Breakdown Structure 
(WBS). 2 
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During the initial stages of developing these technical plans, it was recognized that requirements were needed to 
define "how good" the SLD simulation capabilities needed to be. Therefore, a specific task was incorporated into 
the SLD Engineering Tool Development Project Plan to address this need (Task 1.1). 2 The definition of 
requirements was intended to provide a "target" for SLD simulation in quantified terms. This served the dual 
purpose of providing guidance about: (1) what essential features or characteristics need to be simulated, and (2) how 
accurately these characteristics need to be simulated. It was anticipated that these requirements would be developed 
by means of sensitivity studies, either experimental or computational. 

Two important areas of the project plan, where the definition of requirements would be very beneficial, are: (1) 
updating ice prediction codes to include SLD icing effects, and (2) development and assessment of icing cloud 
instrumentation to measure SLD icing conditions. Though work has been proceeding in these areas, it would be 
desirable to define limits for "how good" the SLD ice prediction capability, and SLD icing condition measurement 
need to be. 

This paper will present the results of an experimental study whose objective was to investigate the feasibility of 
using changes in ice shape features as a potential method for defining requirements. The authors recognized that 
this was a difficult, complex issue to address. Therefore, this first investigation was limited to a very basic 
sensitivity study where changes in LWC, MVD, and temperature were assessed based on perceivable changes in ice 
shape. 

This parameter sensitivity study was conducted at the NASA Glenn Research Center's Icing Research Tunnel 
(IRT) using a 36 inch chord (0.91 m) NACA-0012 airfoil. Liquid Water Content (LWC), drop size (MVD), and 
tunnel static temperature were varied about a nominal value, and the effects of the parameter change on ice shape 
were noted. The experimental method, test setup and test conditions will be described first. Then selected ice 
tracings will be presented and compared on the basis of qualitative and quantitative factors. The report will 
conclude with an assessment of the test method, significance of the results, and recommendations for future 
research. 


II. Experimental Methodology 


A. Approach 

The objective of this experimental sensitivity study was to investigate the feasibility of using changes in ice 
shape features as a means of developing requirements for improvements to SLD ice prediction codes, and the 
instrumentation used to measure SLD icing conditions. It was also desired to assess the practicality of implementing 
such a method, and to document the relative level of ice feature changes corresponding to the LWC, MVD, and 
temperature changes selected for this study. 

Since this was recognized to be a non-trivial task, this sensitivity study was planned to be the first in a sequence 
of experimental investigations. The plan was to use an icing wind tunnel to vary LWC, MVD, and temperature 
conditions about a nominal condition, document changes in the resultant ice shapes, and in some manner use these 
changes as a metric for developing accuracy estimates. Ultimately, we would like to identify the minimum change 
in a parameter (threshold) which produces an observable change in the resulting ice shape. This "threshold" value 
could then be used to develop an estimate of the required accuracy for that parameter. 

Our plan was to keep this first test simple in scope, by choosing an icing spray condition having a nominal LWC, 
MVD, and temperature. We were not certain about what the "threshold" value would be for LWC, MVD, and 
temperature. Therefore, we used our judgment to select changes in LWC, MVD, and temperature which would 
yield an observable change in the resulting ice shapes. This was done with the assumption that the actual 
"threshold" values would probably be less than what was chosen, such that the parameter variations used in this test 
would constitute an upper bound for the actual "threshold" values. 

B. Test Configuration 

This sensitivity study was conducted in the NASA Glenn Icing Research Tunnel (IRT) which is a closed-loop 
refrigerated wind tunnel. The test section is 6 ft (1.83 m) high and 9 ft (2.74 m) wide, and contains a turntable 
assembly which allows for model angle-of-attack changes. A 5000 hp fan provides airspeeds up to 390 mph (empty 
test section). The refrigeration heat exchanger can control air temperature from +40 °F (4.4 °C) to -20 °F (-28.9 °C). 
The water spray system has been calibrated for simulating icing clouds with droplet MVD of 10 to 216pm, and 
Liquid Water Content (LWC) of 0.2 - 3.4 g/m 3 . Figure 1 shows a schematic view of the IRT. 
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Figure 1. NASA Glenn Icing Research Tunnel 


The test article was a 3 foot (0.914 m) chord NACA-0012 
airfoil having a 6 foot (1.83 m) span. It was mounted 
vertically in the test section, and was fixed at a 2.5° angle-of 
attack for all the conditions tested. A heater mat was 
installed on the model to assist in removing the ice accretion 
and making mass measurements. The heater mat extended 
1 1 inches (0.28 m) above and below the vertical centerline of 
the tunnel (which was at 3 feet (0.914 m) from the floor), and 
wrapped around 3.5 inches (0.089 m) on the upper surface of 
the airfoil, and 8.5 inches (0.216 m) on the lower surface. 
The 2.5° AOA was chosen to provide some difference 
between the upper and lower surface ice shapes of this 
symmetrical airfoil. The test article is shown installed in the 
IRT test section in Fig. 2. 

C. Test Conditions 

The parameters of primary interest in this study were 
LWC, MVD, and temperature. Since it was recognized that 
temperature can have a significant effect on the ice shape 
features, depending on the icing regime (e.g. - glaze, mixed, 
rime), it was decided to evaluate the effect of all three 
parameters at two freezing fractions of 0.3, and 0.7. 
Additionally, it was decided to conduct the sensitivity study 
at two airspeeds of 130 knots and 180 knots, for the purpose 
of evaluating the relative significance of airspeed as a 
parameter. Since the variation of LWC, MVD, and 

temperature on the resulting ice shape was evaluated at two 
airspeeds, and two freezing fractions, one could essentially 
organize the sensitivity study into 4 separate sub-tests 



Figure 2. NACA-0012 airfoil installed in the 
test section of the NASA IRT 
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distinguished by freezing fraction and airspeed: r| = 0.3 & V = 130 knots, r| = 0.3 & V = 180 knots, r| = 0.7 & V = 
130 knots, and r| = 0.7 & V = 180 knots. The test conditions for each of the four sub-test categories are listed below: 

Sub-Test #1: 

LWC = 0.6 g/m 3 +/- 0. lg/m 3 , MVD = 80pm +/- 50pm, T s = 19.5°F +/- 3°F for p = 0.3, V= 130 knots 

Sub-Test #2: 

LWC = 0.6 g/m 3 +/- 0.1g/m 3 , MVD = 80pm +/- 50pm, T s = 15.5°F +/- 3°F for p = 0.3, V= 180 knots 

Sub-Test #3: 

LWC = 0.6 g/m 3 +/- 0.1g/m 3 , MVD = 80pm +/- 50pm, T s = 5.5°F +/- 3°F for p = 0.7, V= 130 knots 

Sub-Test #4: 

LWC = 0.6 g/m 3 +/- 0.1g/m 3 , MVD = 80pm +/- 50pm, T s = -0.5°F +/- 3°F for p = 0.7, V = 180 knots 

Supporting rationale for these parameter choices are further described below: 

1. LWC 

Flight research aircraft have measured nominal LWC of 0.3 g/m 3 in SLD icing conditions, with maximum 
values ranging to 0.5 g/m 3 . 3 The nominal value of 0.6 g/m 3 was selected in an attempt to select a nominal LWC as 
close to flight results as possible, yet still remain within the IRT's capability. The variation of +/- 0.1 g/m 3 was 
chosen based on the thought it might be representative of a 15 percent LWC measurement error, which would 
correspond to a practical uncertainty associated with hot-wire LWC instrumentation. 4 5 

2. MVD 

The MVD variation was chosen on a slightly different basis than LWC or temperature variations. It was not 
related to an instrumentation uncertainty, as were LWC and temperature variations, but was selected to allow a wide 
range of MVD variation. A nominal MVD of 80pm was selected to be within the large droplet portion of the IRT 
spray capability. Then a 50pm parameter variation was chosen to allow drop size variation into the Appendix C 
regime, and also to larger drop sizes within the IRT's large droplet capability. Due to facility constraints trying to 
simultaneously match MVD, LWC, and airspeed, it was necessary to select 28pm instead of 30pm, for the minimum 
drop size in this study. It should be noted that the 50pm parameter variation constituted 60 percent of the nominal 
value. 

3. Temperature 

Since it was recognized that temperature had a very significant effect on ice shape features, it was decided to 
conduct the sensitivity study at a freezing fraction r| of 0.3 and 0.7. The concept of freezing fraction, first proposed 
by Messinger, provides a convenient mechanism to account for the thermodynamics of the different icing regimes 
on the leading edge of an airfoil (r| =1.0 .. rime, r)~0.1 .. glaze). For this study, it was decided to select one freezing 
fraction more toward the rime category, and the other freezing fraction more toward the glaze category. It should be 
mentioned that when a parameter was varied, the freezing fraction was not maintained at a value of 0.3, or 0.7. 
Those freezing fractions only applied to the nominal condition. Flowever for convenience, results from variant 
conditions have been associated with the freezing fraction of the nominal condition when listed in tables or referred 
to in figures. 

The nominal values of Ts for this sensitivity study were determined with the aid of LEWICE 2.0 using the 
following procedure: 

• First, the nominal LWC, MVD, airspeed, and AOA and an initial static temperature T s guess (to yield 
the desired freezing fraction) were input to LEWICE 2.0 

• LEWICE was run with those inputs for a 1 minute time period. 

• The value of freezing fraction r| computed for the stagnation point was then compared with the desired 
freezing fraction 

• If the computed freezing fraction differed from the desired value, T s was adjusted and LEWICE was re- 
run again to compute a new freezing fraction at the stagnation point of the test article 

• This process was repeated until the desired and computed freezing fractions agreed 

The 3 °F variation was chosen based on what was believed to be an estimate of a worst case end-to-end 
uncertainty associated with a thermocouple based temperature measurement system (where end-to-end refers to the 
complete signal path from sensor -to- signal conditioning -to- computer output). This estimate assumes +/- 1 °F 
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variation attributed to wiring and signal conditioning, and +/- 2 °F variation for sensor accuracy and artifacts 
introduced in the measurement by the specific sensor installation. 

D. Test Procedure 

The following test procedure was used for each of the test runs in this investigation. 


• Baseline and variation conditions were used to set tunnel operating values (i.e. tunnel velocity, air 
temperature, spray bar air and water pressures) 

• Icing conditions were applied for the selected spray times 

• After the spray was concluded and the tunnel velocity was set to idle: 

• Photographs of the ice accretion were taken 

• An ice shape tracing was made at the mid-span location 

• The ice was cut at predetermined span-wise locations above and below mid-span 

• Heaters were turned on to de-bond the ice from the surface 

• Ice was removed and weighed 

Table 1. Sensitivity study experimental icing runs 


Description 

Run 

n 

V 

(knots) 

LWC 

(g/m 3 ) 

MVD 

(pm) 

Ttot 

(°F) 

Ts 

(°F) 

AOA 

(degrees) 

Mass 

(gm) 

Spray 

Time 

(min) 

Pair 

(psi) 

AP 

(psi) 

nominal 

1 

0.3 

130 

0.6 

80 

23.5 

19.5 

2.5 

845.6 

20 

4.5 

19.1 

repeat run 

2 

0.3 

130 

0.6 

80 

23.5 

19.5 

2.5 

882.9 

20 

4.5 

19.1 

LWC - 0.1 

3 



0.5 

" 

" 

" 

t! 

740.9 

II 

3.4 

12.4 

LWC + 0.1 

4 



0.7 

" 

" 

" 

t! 

1027.5 

II 

6.0 

28.2 

MVD - 52 

5 



0. 6 

28 

" 

" 

It 

715.9 

II 

10.1 

30.3 

MVD + 50 

6 



" 

130 

" 

" 

II 

1081.2 

II 

2.5 

14.8 

Temp - 3 

7 



" 

80 

20.5 

16.5 

II 

862.2 

II 

4.5 

19.1 

Temp + 3 

9 



" 

" 

26.5 

22. 5 

II 

877.9 

II 

4.5 

19.1 

nominal 

10 

0.7 

130 

0.6 

80 

9.5 

5.5 

2.5 

1059.9 

20 

4.5 

19.1 

repeat run 

11 

0.7 

130 

0.6 

80 

9.5 

5.5 

2.5 

1102.5 

20 

4.5 

19.1 

LWC-0.1 

12 



0.5 

" 

" 

" 

II 

970.7 

II 

3.4 

12.4 

LWC + 0.1 

13 



0.7 

" 

" 

" 

II 

1241.9 

II 

6.0 

28.2 

MVD - 52 

14 



0. 6 

28 

" 

" 

II 

838.9 

II 

10.1 

30.3 

MVD + 50 

15 



" 

130 

" 

" 

II 

1458.8 

II 

2.5 

14.8 

Temp - 3 

16 



" 

80 

6.5 

2. 5 

II 

1168.2 

II 

4.5 

19.1 

Temp + 3 

17 



" 

" 

12.5 

8. 5 

II 

1058.1 

II 

4.5 

19.1 

nominal 

18 

0.3 

180 

0.6 

80 

23.2 

15.5 

2.5 

784.2 

13.6 

8.2 

42.2 

repeat run 

19 

0.3 

180 

0.6 

80 

23.2 

15.5 

2.5 

763.5 

13.6 

8.2 

42.2 

LWC-0.1 

20 



0.5 

" 

" 

" 

II 

649.6 

II 

5.6 

25.8 

LWC + 0.1 

21 



0.7 

" 

" 

" 

it 

1016.5 

II 

10 

65 

MVD - 52 

22 



0. 6 

28 

" 

" 

II 

611.2 

II 

13.2 

47.9 

MVD + 50 

23 



" 

130 

" 

" 

II 

844.8 

II 

3.7 

27 

Temp - 3 

24 



" 

80 

20.2 

12.5 

II 

785.1 

II 

8.2 

42 

Temp + 3 

25 



" 

" 

26.2 

18.5 

II 

755.1 

II 

8.2 

42 

nominal 

27 

0.7 

180 

0.6 

80 

7.2 

-0.5 

2.5 

1009.1 

13.6 

8.2 

42.2 

repeat run 

28 

0.7 

180 

0.6 

80 

7.2 

-0.5 

2.5 

1014.3 

13.6 

8.2 

42.2 

LWC-0.1 

29 



0.5 

" 

" 

" 

II 

818.9 

II 

5.6 

25.8 

LWC + 0.1 

30 



0.7 

" 

" 

" 

II 

1424.3 

II 

10 

65 

MVD - 52 

31 



0. 6 

28 

" 

" 

II 

684.1 

II 

13.2 

47.9 

MVD + 50 

32 



" 

130 

" 

" 

it 

1275.4 

II 

3.7 

27 

Temp - 3 

33 



" 

80 

4.2 

-3.5 

II 

1099.1 

II 

8.2 

42 

Temp + 3 

34 



" 

" 

10.2 

2.5 

II 

1058.6 

II 

8.2 

42 


Note: accumulation parameter was held constant for the two airspeeds, which resulted in different spray times 


NASA/TM— 2005-213562 


5 


TIT. Results and Discussion 

The experimental icing runs conducted in the NASA-Glenn IRT as part of this sensitivity investigation are listed 
in Table 1. The tunnel and icing spray conditions are tabulated for each run, along with the mass of the ice accretion 
removed from the area of the model covered by the heater mat. This section will present selected results, and 
describe the qualitative and quantitative analysis applied to these results. The ice shape repeatability will be 
considered, followed by an assessment of LWC, MVD, and temperature variations on the resulting ice shapes. 
Though the study was done at two airspeeds, space limitations allow us to present only the data for one airspeed. 
Since, the higher airspeed was closer to actual aircraft speeds, it was thought these results might have more 
relevance to aircraft icing. Therefore, only the 180 knot cases will be presented and discussed. 

A. Ice Shape Repeatability 

There is an inherent level of variation to ice shapes produced in an icing wind tunnel. Since this study intended 
to evaluate the effects of parameter changes on the resulting ice accretion, it was important to determine what level 
of change in ice shape features (general shape, horn location, mass... etc) could be expected due to the normal run to 
run variations of the icing spray condition. Obviously, less ice shape variation (better repeatability) would be more 
desirable because it would facilitate the evaluation of the effects of smaller parameter changes. 



X, in 

c) n = 0.3, V = 130 knots, Ts = 19.5 °F , 
LWOO.6 g/m 3 , MVD = 80pm, AOA = 2.5° 



X, in 

d) ii = 0.3, V = 180 knots, Ts = 15.5 °F, 

LWC = 0.6 g/m 3 , MVD = 80pm, AOA = 2.5° 



X, in 

a) q = 0.7, V = 130 knots, Ts = 5.5 °F 

LWC = 0.6 g/m 3 , MVD = 80pm, AOA = 2.5° 



X, in 

b) q = 0.7, V = 180 knots, Ts = -0.5 °F 

LWC = 0.6 g/m 3 , MVD = 80pm, AOA = 2.5° 


Figure 3. Ice shape repeatability for repeat runs of the 4 nominal sub-test conditions. 
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1. Qualitative Evaluation 

Figure 3 shows ice tracing comparisons from repeat runs of the nominal spray condition having an LWC of 
0.6 g/m 3 and an MVD of 80pm. These ice tracings were obtained at two freezing fractions (r| = 0.3 and 0.7), and two 
airspeeds (V=130 and 180 knots). A qualitative inspection of the tracing comparison plots in fig. 3 suggests that the 
repeatability was very good. Even at warmer temperature conditions (r| = 0.3), there was very good agreement 
between the main ice shapes for both airspeed conditions (fig. 3c and 3d). The ice tracing comparison in fig. 3a 
shows more noticeable variation than figs. 3b, 3c, and 3d, but this is not considered significant because there may 
have been an issue with the ice tracing for Run 10 (inner tracing). It is believed that the cardboard ice tracing 
template was not pushed all the way into the ice cut on the accretion for this run, resulting in the observed variation. 


2. Quantitative Evaluation 

It was also desired to employ some quantitative means to assess changes in ice shape features in addition to 
qualitative comparisons. The THICK program developed by William Wright of NASA Glenn Research Center was 
used to perform quantitative assessment of changes in ice shape features. 6 This program was originally developed 
and used during the validation of LEWICE. By comparing the un-iced airfoil coordinates with the digitized 
coordinates of an ice shape, it is capable of calculating ice shape attributes such as upper and lower surface 
maximum ice thickness, icing horn angle, 
and icing limits. 

Fig. 4 shows the convention that 
THICK uses to determine maximum 
thickness and horn angle. An axis has been 
superimposed on the ice tracing to illustrate 
the coordinate system used by THICK to 
determine the angles of the maximum ice 
thickness (which THICK considers as 
horns) on both the upper and lower surface 
of the airfoil. The maximum thickness is 
calculated by scanning the ice shape for the 
location of greatest thickness normal to the 
surface on both the upper and lower 
surface of the airfoil. When a maximum 
thickness is found, THICK then draws a 
line from the origin of the axis system to 
the point of maximum thickness, 
whereupon it then calculates the associated 
angle as depicted in Fig. 4. 

THICK was used to determine the ice horn thickness, horn angle, iced area, and icing limits associated with the 
ice tracings shown in Figs. 3a-3d. Since each plot in Fig. 3 contained two ice tracings (representing repeat runs of 

Table 2. THICK Calculations For Repeatability Of Nominal Sub-Test Conditions 



x (in) 

Figure 4. Diagram showing THICK computation of 


Run 

Freezing 

Fraction V MVD 

H (knots) (pm) 

LWC, 

(g/m 3) 

Lower 
Icing Limit 
(inches) 

Upper Icing 
Limit 
(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass (gm) 

1 

0.3 130 80 

0.6 

-7.5 

2.36 

1.33 

1.61 

4.34 

285.9 

120.8 

845.6 

2 



-7.77 

3.14 

1.25 

1.50 

4.27 

281.6 

113.5 

882.9 

AVG 1&2 

Nominal Sub-Test Condition #1 

-7.66 

2.75 

1.29 

1.56 

4.31 

283.8 

117.2 

864.3 

STD DEV 

Variation Of Sub-Test Condition #1 

0.16 

0.55 

0.06 

0.08 

0.05 

3.0 

5.2 

26.4 

18 

0.3 180 80 

0.6 

-5.87 

1.83 

1.40 

2 

4.20 

288.5 

126.1 

784.2 

19 



-6.30 

2.00 

1.56 

1.44 

3.97 

286.3 

123.9 

763.5 

AVG 18&19 

Nominal Sub-Test Condition #2 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 

STD DEV 

Variation Of Sub-Test Condition #2 

0.31 

0.12 

0.12 

0.06 

0.17 

1.6 

1.6 

14.6 

10 

0.7 130 80 

0.6 

-7.70 

3.34 

1.19 

1.21 

5.33 

236.7 

124.4 

1059.9 

11 



-6.43 

3.09 

1.56 

1.52 

5.70 

230.9 

137.8 

1102.5 

AVG 10 & 11 

Nominal Sub-Test Condition #3 

-7.07 

3.22 

1.37 

1.37 

5.52 

233.8 

131.1 

1081.2 

STD DEV 

Variation Of Sub-Test Condition #3 

0.90 

0.18 

0.26 

0.22 

0.26 

4.1 

9.5 

30.1 

27 

0.7 130 80 

0.6 

-7.67 

3.03 

1.63 

1.78 

5.72 

225.5 

172.1 

1009.1 

28 



-7.70 

3.06 

1.48 

1.59 

5.51 

225.2 

172.5 

1014.3 

AVG 27 & 28 

Nominal Sub-Test Condition #4 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 

STD DEV 

Variation Of Sub-Test Condition #4 

0.03 

0.02 

0.11 

0.13 

0.15 

0.2 

0.3 

3.7 
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the same nominal sub-test condition), THICK was run twice to get a unique value of the previously mentioned ice 
shape attributes. The measured mass was also included with the corresponding THICK calculations. These values 
were then used to determine an average and standard deviation of each ice shape attribute for each nominal sub-test 
condition. These results are listed in Table 2, for each of the nominal sub-test conditions. 

The values in Table 2 will provide the basis for quantitatively assessing the significance of ice feature changes 
due to variations in LWC, MVD, and temperature later in this report. The significance of the average value of an 
attribute in Table 2, is that they will serve as reference or (nominal) THICK values for assessing the effect of LWC, 
MVD, and temperature effects. The standard deviation of an attribute in Table 2 represents the expected variation 
in that ice shape attribute due to icing spray repeatability for the particular nominal sub-test condition. 

Inspection of the standard deviation values in Table 2, suggests that for all THICK attributes (except the icing 
limits) there was a relatively small variation. This is a quantitative indication of good ice shape repeatability for all 
the cases in Fig. 3. William Wright obtained similar results during the validation of LEWICE 2.0. 6 The greater 
variation in the standard deviation of THICK icing limits is not totally unexpected, given the difficulties of tracing 
ice features near the icing limits. It is sometimes very difficult to trace small feathers far aft of the main ice shape, 
without them breaking off. 

B. Effect of Varying LWC 



X, in 

a) p = 0.3, V = 180 knots, Ts = 15.5 °F, 
MVD = 80pm, AOA = 2.5° 



X, in 

b) p = 0.7, V = 180 knots, Ts = -0.5 °F 
MVD = 80pm, AOA = 2.5° 


Figure 5. LEWICE predictions showing the effect of LWC variations on ice shape. 



X, in 


a) p = 0.3, V = 180 knots, Ts = 15.5 °F, 
MVD = 80pm, AOA = 2.5° 



X, in 

b) p = 0.7, V = 180 knots, Ts = -0.5 °F 
MVD = 80pm, AOA = 2.5° 


Figure 6. Change in ice shape tracings due to LWC variation at p = 0.3 and 0.7 
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1. Computational Trend Prediction 

Since, it was known that freezing fraction had a significant effect on ice shape features, the effect of LWC 
variations on ice shapes was evaluated at two freezing fractions of 0.3 and 0.7. To get a sense of how ice shape 
features might change as a result of LWC variations, LEWICE 2.0 was run with the appropriate input conditions 
from Table 1. The predicted ice shapes for a 180 knot airspeed are shown in Fig. 5. 

Figure 5 illustrates the effect of LWC variation on predicted ice shapes for LWC variations of +/- 0.1 g/m 3 about 
a nominal condition of 0.6 g/m 3 . The effect of freezing fraction on ice shape feature changes can be seen by 
comparing Fig. 5a where r| = 0.3, with Fig. 5b where r| = 0.7. Some qualitative observations can be made about 
expected ice shape changes based on Fig. 5. 

For a freezing fraction of 0.3, it appears as if: 

• The thickness of the flat area main ice shape remains the same regardless of LWC. 

• As LWC is increased the upper horns grow larger, and rotate aft (decrease in TFIICK horn angle). 

• The main ice shape extends further aft on the lower surface as the LWC is increased. 

• The mass and area of the main ice shape increase as LWC is increased. 

For a freezing fraction of 0.7, it appears as if: 

• Ice horn thickness increases as LWC increases. 

• The upper ice horn angle doesn't appear to change much. 

• The lower ice horn angle appears to rotate aft a bit. 

• The area of the main ice shape increases with increasing LWC. 

The LEWICE results suggest that the +/- 0.1 g/m 3 LWC variation should yield discernible changes in ice shape 
features, and that the type and degree of feature changes appear to be affected by the freezing fraction. 

2. Qualitative Assessment Of Actual Ice Shapes 

Fig. 6 illustrates the effect of LWC variation on actual ice shapes for LWC variations of +/- 0.1 g/m 3 about a 
nominal condition of 0.6 g/m 3 . The effect of freezing fraction on ice shape feature changes can be seen by 
comparing Fig. 6a where r| = 0.3, with Fig. 6b where r| = 0.7. 

Examination of the ice shape tracings in Fig. 6a and Fig. 6b exhibited similar trends to what was predicted by 
LEWICE. Since LEWICE doesn't model feather growth, there were some expected differences between the actual 
and predicted ice shapes. Some of these differences between predicted and observations about actual ice shape 
changes are listed below. 

For a freezing fraction of 0.3: 

• The thickness of the flat area of the main ice shape remained about constant (regardless of LWC) near 
the stagnation line, but the thickness began to change as a function of LWC as you moved toward the 
lower surface 

• As LWC was increased, the lower surface horns became thicker, and appeared to rotate aft. The 
differences in horn thickness and rotation appeared to be greater when comparing the nominal condition 
(0.6 g/m 3 ) to the high LWC condition (0.5 g/m 3 ), versus comparing the nominal condition to the low 
LWC condition (0.5 g/m 3 ). 

• The lower icing limit appeared to move forward as LWC was increased 
For a freezing fraction of 0.7: 

• In general, there seemed to be less relative change in ice shape features as LWC was varied, than the 
LEWICE predictions. 

• The presence of large ice feathers was evident on the actual ice tracings at this freezing fraction. 

• Both the upper and lower ice horn angles didn't appear to change significantly with changes in LWC. 
Based on the above qualitative assessment, it can be said that a +/- 0.1 g/m 3 variation in LWC produced visually 
observable changes in ice shape features, and that the freezing fraction affected the degree to which these feature 
changes were manifested. 

3. Quantitative Assessment Of Actual Ice Shapes 

Given the qualitative assessment above, it was next desired to quantitatively evaluate changes in ice shape 
features due to variations in LWC. The TFIICK program was used to quantify ice shape features and the change in 
these features as the LWC was varied +/- 0.1 gm/m 3 about the nominal LWC of 0.6 gm/m 3 . Tables 3 and 4 contain 
the results of the THICK analysis for r| = 0.3 and r| = 0.7 respectively. Each table shows the THICK calculated ice 
shape attributes for the low LWC condition 0.5 g/m 3 , the nominal LWC condition 0.6 g/m 3 , and the high LWC 
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condition 0.7 g/m 3 on the top 3 rows of the table. The ice shape tracings in Figs. 5 and 6 were input to the THICK 
program to generate these values. Though the mass was included in these tables, it was directly measured by 
removing the ice accretion over a known part of the airfoil, and then directly weighing it. Thus it was not calculated 
by THICK. 

It should be noted that the THICK calculations for the nominal condition are actually an average of THICK 
values from repeat runs of the nominal condition (this was discussed earlier in the section of this paper dealing with 
repeatability). In addition, a standard deviation was determined from the THICK calculations for those repeat runs. 

The standard deviation was intended to provide a "threshold" value by which to judge if the change in a THICK 
ice shape attribute could be considered significant. The rationale for this was that the standard deviation of an ice 
shape attribute essentially represented it's inherent variability due to variations in the icing spray condition from run 
to run (i.e., spray condition repeatability). The standard deviation for each ice shape attribute is listed in the second 
row from the bottom of each table. 

The lower 3 rows of each table show the calculated change of an ice shape attribute with respect to it's value at 
the nominal condition. The absolute value of these changes was compared with the standard deviation, and 
wherever the change exceeded the standard deviation, it was highlighted. It should be mentioned that the following 
convention applies to the data in the tables below: (1) positive change in icing limit on upper surface = aft 
movement, (2) positive change in icing limit on lower surface = forward movement, and (3) + change in horn angle 
= counter-clockwise rotation. 

Table 3. Thick Calculations For Ice Shape Changes Due To LWC Variation At ?i = 0.3 


Description 

Run 

n 

V 

(knots) 

LWC, 

(g/m 3) 

Lower 

Icing 

Limit 

(inches) 

Upper 

Icing 

Limit 

(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass 

(gm) 

THICK 

20 



0.5 

-7.5 

2.029 

1.340 

1.248 

3.65 

279.7 

131.6 

649.6 

Values 

Nominal 

0.3 

180 

0.6 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 


21 



0.7 

-5.23 

1.97 

1.49 

1.620 

4.79 

298.9 

119.6 

1016.5 

Delta 

20 



-0.1 

- 1.42 

0.11 

- 0.14 

- 0.232 

- 0.44 

- 7.7 

6.6 

- 124.3 

WRT 

Nominal 

0.3 

180 

0.0 

0.00 

0.00 

0.00 

0 

0.00 

0.0 

0.0 

0.0 

Nominal 

21 



0.1 

0.86 

0.06 

0.01 

0.14 

0.71 

11.5 

- 5.4 

242.7 



Nominal 

Runs 18&19 



Std Dev 

0.31 

0.12 

0.12 

0.06 

0.17 

1.6 

1.6 

14.6 

Condition 

Runs 18&19 



Average 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 


Table 4. Thick Calculations For Ice Shape Changes Due To LWC Variation At = 0.7 


Description 

Run 

n 

V 

(knots) 

LWC, 

(g/m 3) 

Lower 

Icing 

Limit 

(inches) 

Upper 

Icing 

Limit 

(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass 

( gm ) 

THICK 

29 



0.5 

-7.73 

1.94 

1.38 

1.40 

4.81 

228.0 

178.9 

818.9 

Values 

Nominal 

0.7 

180 

0.6 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 


30 



0.7 

-7.67 

2.16 

1.67 

1.87 

6.56 

220.2 

171.0 

1424.3 

Delta 

29 



-0.1 

- 0.04 

- 1.11 

- 0.17 

- 0.29 

- 0.81 

2.7 

6.6 

- 192.8 

WRT 

Nominal 

0.7 

180 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

0.0 

0.0 

Nominal 

30 



0.1 

0.02 

- 0.89 

0.12 

0.18 

0.94 

- 5.1 

- 1.3 

412.6 



Nominal 

Runs 27&28 



Std Dev 

0.03 

0.02 

0.11 

0.13 

0.15 

0.2 

0.3 

3.7 

Condition 

Runs 27&28 



Average 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 


THICK values in Table 3 were compared with the ice tracings in Fig. 6a. This comparison revealed that except 
for the mass measurements, the quantitative changes tracked the qualitative observations quite well. Likewise 
THICK values from Table 4 were also compared with the ice tracings in Fig. 6b, and similar results were found. 

For some unexplained reason the change in mass was greater for the LWC increase than for the LWC decrease. 
A review of the mass change from similar conditions at 130 knots, showed that the mass changes were 
approximately the same for an increase and decrease in LWC. Since the mass change associated with the LWC 
increase in Tables 3 and 4 appeared to be rather large, we suspected there might be a problem with the reported 
LWC for that condition. Therefore, we plan to run an icing blade at the condition in question to verify the actual 
LWC. 
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There was also an apparent issue with the THICK calculated value of lower horn thickness for run 29 in Table 4. 
This value of -0.1 7 agreed with the general trend for the lower horn thickness to decrease with decreasing LWC. 
However, the actual value does not compare well with the observed change in Fig. 6b, where the lower horn peak 
for run 29 just about overlays the peak for run 28 (which was one of the nominal condition runs). Based on this 
observation the change in lower ice horn thickness would be negligible. 

The reason for the apparent discrepancy between the observed change in lower horn thickness in Fig. 6b, and the 
calculated value of -0.0 17 in Table 4, is that nominal values in Table 4 were calculated using the average value of 
lower horn thickness from runs 27 and 28. If we had plotted an ice shape that represented the average of runs 27 
and run 28, the calculated change would agree. 

Though this appears to be a rather insignificant issue, this apparent discrepancy illustrates the importance of 
having an accurate nominal or reference value to compare quantitative ice feature changes against. Therefore, it 
would be desirable to use THICK to determine the nominal ice shape attribute values from more than two repeat 
runs. 

C. Effect of Varying MVD 




a) ii = 0.3, V = 180 knots, Ts = 15.5 °F, 
LWC = 0.6 g/m 3 , AOA = 2.5° 


b) n = 0.7, V = 180 knots, Ts = -0.5 °F 
LWC = 0.6 g/m 3 , AOA = 2.5° 


Figure 7. LEWICE predictions showing the effect of MVD variation on ice shape. 




X, in X, in 


a) n = 0.3, V = 180 knots, Ts = 15.5 °F, 
LWC = 0.6 g/m 3 , AOA = 2.5° 


b) n = 0.7, V = 180 knots, Ts = -0.5 °F 
LWC = 0.6 g/m 3 , AOA = 2.5° 


Figure 8. Change in ice shape tracings due to MVD variation at = 0.3 and 0.7 
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1. Computational Trend Prediction 

The effect of MVD variations on ice shape features was evaluated at two freezing fractions of 0.3 & 0.7. To get a 
sense of how ice shape features might change as a result of MVD variations, LEWICE 2.0 was run with the 
appropriate input conditions from Table 1. The predicted ice shapes for a 180 knot airspeed are shown in Fig. 7. 

Figure 7 illustrates the effect of MVD variation on predicted ice shapes for MVD variations of +/- 50pm about a 
nominal condition of 80pm. The effect of freezing fraction on ice shape feature changes can be seen by comparing 
Fig. 7a where r| = 0.3, with Fig. 7b where r| = 0.7. Some qualitative observations can be made about expected ice 
shape changes based on Fig. 7. 

For a freezing fraction of 0.3, it appears as if: 

• There is a relatively flat area of the main ice shape about the stagnation region which remains 
approximately the same regardless of MVD. 

• As MVD is increased the upper horns grow larger, and rotate aft (decrease in TF1ICK horn angle). 

• The main ice shape extends further aft on the lower surface as the MVD is increased. 

• The mass and area of the main ice shape increase as MVD is increased. 

• Generally speaking, there is a much smaller change in the ice shape features for an increase in MVD to 
130pm relative to the nominal test condition of 80pm, than for a decrease in MVD to 28pm. 

For a freezing fraction of 0.7, it appears as if: 

• Ice horn thickness increases as MVD increases. 

• Ice horn angles don't appear to change much. 

• There is more ice aft of the main ice shape on both surfaces as MVD is increased. 

• As with the r| = 0.3 observations, there is a much smaller change in the ice shape features for an 
increase in MVD to 130pm relative to the nominal test condition of 80pm, than for a decrease in MVD 
to 28 pm. 

2. Qualitative Assessment Of Actual Ice Shapes 

Figure 8 illustrates the effect of MVD variation on actual ice shapes for MVD variations of +/- 50pm about a 
nominal condition of 80pm. The effect of freezing fraction on ice shape feature changes can be seen by comparing 
Fig. 8a where r| = 0.3, with Fig. 8b where r| = 0.7. 

The ice shape tracings in Figs. 8a and 8b exhibited similar trends to what was predicted by LEWICE. However, 
there were some differences between the actual and predicted ice shapes, and these are listed below. 

For a freezing fraction of 0.3: 

• The thickness of the flat area of the main ice shape remained about constant (regardless of MVD) near 
the stagnation line, but the thickness began to change as a function of MVD as you moved toward the 
lower surface. 

• The appearance of large glaze feathers became more prominent aft of the main ice shape as MVD was 
increased 

• As MVD was increased, the lower surface horns became thicker, and appeared to rotate aft. Unlike the 
LEWICE predictions, there seemed to be an equally significant change in lower horn thickness for a 
decrease in MVD, when comparing the change observed for an increase in MVD. 

• The upper icing limit appeared to move afl on the upper surface for the 130pm run. 

• The lower icing limit appeared to move forward on the lower surface with increased MVD 

For a freezing fraction of 0.7: 

• In general, the observed trends seemed to follow the LEWICE prediction trends. 

• As MVD was increased the feathers on the main ice shape became more clearly defined. 

• The icing limit for the 130pm condition extended further aft on the upper surface than for the 28pm, or 
80pm 

Based on the above qualitative assessment, it can be said that a +/- 50pm variation in MVD produced visually 
observable changes in ice shape features, and that the freezing fraction affected the degree to which these feature 
changes were manifested. 

3. Quantitative Assessment of Actual Ice Shapes 

Given the qualitative assessment above, it was next desired to quantitatively evaluate changes in ice shape 
features due to variations in MVD. The THICK program was used to quantify ice shape features and the change in 
these features for an MVD variation of +/- 50pm. Tables 5 and 6 contain the results of the THICK analysis for r| = 
0.3 and r| = 0.7, respectively. 
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Table 5. Thick Calculations For Ice Shape Changes Due To MVD Variation At ij = 0.3 


Description 

Run 

n 

V 

(knots) 

MVD, 

(pm) 

Lower 

Icing 

Limit 

(inches) 

Upper 

Icing 

Limit 

(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass 

(gm) 

THICK 

22 



28 

-7.63 

1.92 

1.35 

1.46 

3.32 

272.0 

141.0 

611.2 

Values 

Nominal 

0.3 

180 

80 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 


23 



130 

-5.42 

3.01 

1.26 

1.72 

5.18 

299.2 

121.2 

844.8 

Delta 

22 



-52 

- 1.55 

0.01 

- 0.13 

-0.02 

- 0.76 

- 15.4 

16.0 

- 162.7 

WRT 

Nominal 

0.3 

180 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

0.0 

0.0 

Nominal 

23 



50 

0.67 

1.10 

- 0.22 

0.24 

1.09 

11.8 

- 3.8 

70.9 



Nominal 

Runs 18&19 



Std Dev 

0.31 

0.12 

0.12 

0.06 

0.17 

1.6 

1.6 

14.6 

Condition 

Runs 18&19 



Average 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 


Table 6 Thick Calculations For Ice Shape Changes Due To MVD Variation At r| = 0.7 


Description 

Run 

n 

V 

(knots) 

MVD, 

(pm) 

Lower 

Icing 

Limit 

(inches) 

Upper 

Icing 

Limit 

(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass 

(gm) 

THICK 

31 



28 

-7.67 

1.92 

1.16 

1.25 

3.96 

215.1 

186.6 

684.1 

Values 

Nominal 

0.7 

180 

80 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 


32 



130 

-7.58 

1.95 

1.82 

1.95 

7.78 

221.4 

169.7 

1275.4 

Delta 

31 



-52 

0.02 

- 1.13 

- 0.39 

- 0.43 

- 1.65 

- 10.2 

14.3 

- 327.6 

WRT 

Nominal 

0.7 

180 

0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

0.0 

0.0 

Nominal 

32 



50 

0.11 

- 1.10 

0.27 

0.26 

2.17 

- 3.9 

- 2.6 

263.7 



Nominal 

Runs 27&28 



Std Dev 

0.03 

0.02 

0.11 

0.13 

0.15 

0.2 

0.3 

3.7 

Condition 

Runs 27&28 



Average 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 


THICK values from Table 5 were compared with the ice tracings in Fig. 8a. Ice shape feature changes which 
exceeded the standard deviation have been highlighted. This comparison revealed that the quantitative changes 
tracked the qualitative observations quite well. 

There were a few attribute changes that looked questionable, however. These were: (1) the upper icing limit for 
run 23, and (2) the apparent forward movement of the lower icing limits with increased drop size. After reviewing 
photographs, the run 23 upper icing limit change was due to a small surface feature that accreted on a surface 
irregularity, and probably should not have been traced. Thus it can be considered to be an ice tracing artifact, and 
not real. 



a) MVD = 28pm, run 22 b) MVD = 80pm, run 19 c) MVD = 130pm, Run 23 

Figure 9. Comparison of lower surface icing limits for different MVD (note: LE is toward left side of picture). 
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Interestingly enough, the apparent forward movement of the lower surface icing limit was actually real. Figure 9 
illustrates this trend with 3 photos of the lower surface. 

THICK values from Table 6 were compared with the ice tracings in Fig. 8b, and it was found those quantitative 
changes again tracked the qualitative observations. The calculated value for upper icing limit change had a negative 
value for the low MVD condition (run 31) and also for the high MVD condition (run 32). This looked strange, so 
photos of the upper surface were reviewed for these test runs, and it was determined that this was an ice tracing 
artifact, and thus not a real change. 


D. Effect of Varying Temperature 



X, in 



X, in 


a) ii = 0.3, V = 180 knots, LWC = 0.6 gm/m 3 , 
MVD = 80pm, AOA = 2.5° 


b) i\ = 0.7, V = 180 knots, LWC = 0.6 gm/m 3 , 
MVD = 80pm, AOA = 2.5° 


Figure 10. LEWICE predictions showing the effect of temperature variations on ice shape. 



X, in 

a) p = 0.3, V = 180 knots, LWC = 0.6 gm/m 3 , 
MVD = 80pm, AOA = 2.5° 

Figure 11. Change in ice shape tracings 



X, in 

b) r\ = 0.7, V = 180 knots, LWC = 0.6 gm/m 3 , 
MVD = 80pm, AOA = 2.5° 

to temperature variation at p = 0.3 and 0.7 
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1. Computational Trend Prediction 

The effect of temperature variations on ice shapes was evaluated at two freezing fractions of 0.3 and 0.7. To get 
a sense of how ice shape features might change as a result of temperature variations, LEWICE 2.0 was run with the 
appropriate input conditions from Table 1. The predicted ice shapes for a 180 knot airspeed are shown in Fig. 10. 

Figure 10 illustrates the effect of MVD variation on predicted ice shapes for MVD variations of +/- 3 °F about 
nominal conditions of 15.5 °F at r| = 0.3, and -0.5 °F at q = 0.7. The effect of freezing fraction on ice shape feature 
changes can be seen by comparing Fig. 10a where r| = 0.3, with Fig. 10b where r| = 0.7. Some qualitative 
observations can be made about expected ice shape changes based on Fig. 10. 

For a freezing fraction of 0.3, it appears as if: 

• The upper ice horns rotate aft toward the upper surface as the temperature is increased 

• The mass in the main ice shape is re-distributed further aft on the lower surface, as the temperature 
increases. 

For a freezing fraction of 0.7, it appears as if: 

• The upper ice horn thickness and angles remain the same with increasing temperature. 

• The lower ice horn angle rotates aft toward the lower surface with increasing temperature. 

The LEWICE results suggest that the +/- 3 °F variation should yield discernible changes in ice shape features, 
and that the feature changes appear to be affected by the freezing fraction. 

2. Qualitative Assessment Of Actual Ice Shapes 

Figure 1 1 illustrates the effect of temperature variation on actual ice shapes for temperature variations of +/- 3 °F 
about nominal conditions of 15.5 °F at r| = 0.3, and -0.5 °F at r| = 0.7. The effect of freezing fraction on ice shape 
feature changes can be seen by comparing Fig. 11a where r| = 0.3, with Fig. 1 lb where q = 0.7. 

The ice shape tracings in Figs. 11a and lib exhibited similar trends to what was predicted by LEWICE. 
However, there were some differences between the actual and predicted ice shapes, and these are listed below. 

For a freezing fraction of 0.3: 

• The upper ice horns showed a very significant rotation aft toward the upper surface (decreasing horn 
angle). 

• The lower horn thickness for the warmest temperature case seems to have decreased significantly, with 
the re-distribution of mass on the lower surface of the airfoil. 

• The icing limits moved very slightly aft on both the upper and lower surfaces 

For a freezing fraction of 0.7: 

• In general, the horn thickness appears to have decreased with increasing temperature 

• The upper horn angle appears to have rotated slightly aft toward the upper surface with increasing 
temperature. 

• The icing limits appear to have moved slightly forward toward the leading edge on both surfaces. 

Based on the above qualitative assessment, it can be said that a +/- 3 °F variation in temperature produced 

visually observable changes in ice shape features, and that the freezing fraction affected the degree to which these 
feature changes were manifested. 

3. Quantitative Assessment of Actual Ice Shapes 

Given the qualitative assessment above, it was next desired to quantitatively evaluate changes in ice shape 
features due to these variations in temperature. The THICK program was used to quantify ice shape features and the 
change in these features for a +/- 3 °F temperature change. Tables 7 and 8 contain the results of the THICK 
analysis for q = 0.3 and q = 0.7, respectively. 

THICK values from Table 7 were compared with the ice tracings in Fig. 11a. Ice shape feature changes which 
exceeded the standard deviation have been highlighted. This comparison revealed that quantitative changes agreed 
well with qualitative observations about ice feature change due to temperature variation. 

There were however several ice shape feature changes in Table 7 that are worthy of mention. One of those was 
the negligible mass change with temperature variation. Since the icing spray condition remained the same for each 
test run in Table 7, it was somewhat expected that there would be a small variation in mass. The mass change 
results in Table 7 were only slightly larger than the standard deviation associated with the nominal condition. This 
result also implies that if the mass isn't changing with temperature, then any ice feature changes would be due to a 
re-distribution of mass. This rationale agrees with the qualitative observations of ice feature changes noted in 
Fig. 11a, which indicated the main ice shape mass seemed to shift progressively more toward the lower surface with 
increasing temperature. 
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It is interesting to note that the ice horn angle changes in Table 7, were greater than any horn angle changes due 
to LWC or MVD variations discussed so far. They were about twice the maximum value noted for ice horn angle 
changes due to LWC or MVD variations. 

Table 7. Thick Calculations For Ice Shape Changes Due To Temperature Variation At r| = 0.3 


Description 

Run 

n 

V 

(knots) 

Ts, (°F) 

Lower 

Icing 

Limit 

(inches) 

Upper 

Icing 

Limit 

(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass 

( gm ) 

THICK 

24 



12.5 

-6.18 

1.93 

1.58 

1.69 

3.82 

265.1 

149.9 

785.1 

Values 

Nominal 

0.3 

180 

15.5 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 


25 



18.5 

-6.61 

2.15 

0.83 

1.31 

3.69 

318.8 

96.4 

755.1 

Delta 

24 



-3.0 

-0.09 

0.02 

0.10 

0.21 

- 0.26 

- 22.3 

24.9 

11.3 

WRT 

Nominal 

0.3 

180 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

0.0 

0.0 

Nominal 

25 



3.0 

- 0.53 

0.23 

- 0.65 

- 0.17 

- 0.39 

31.4 

- 28.6 

- 18.8 



Nominal 

Runs 18&19 



Std Dev 

0.31 

0.12 

0.12 

0.06 

0.17 

1.6 

1.6 

14.6 

Condition 

Runs 18&19 



Average 

-6.09 

1.91 

1.48 

1.48 

4.08 

287.4 

125.0 

773.9 


Table 8. Thick Calculations For Ice Shape Changes Due To Temperature Variation At = 0.7 


Description 

Run 

n 

V 

(knots) 

Ts, (°F) 

Lower 

Icing 

Limit 

(inches) 

Upper 

Icing 

Limit 

(inches) 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Area 

(inch 2 ) 

Lower 

Horn 

Angle 

(deg) 

Upper 

Horn 

Angle 

(deg) 

Mass 

(gm) 

THICK 

33 



-3.5 

-7.68 

3.07 

1.66 

1.72 

5.51 

222.5 

178.9 

1099.1 

Values 

Nominal 

0.7 

180 

-0.5 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 


34 



2.5 

-6.80 

2.20 

1.39 

1.49 

4.41 

224.0 

161.6 

1058.6 

Delta 

33 



-3.0 

0.01 

0.02 

0.10 

0.03 

-0.11 

- 2.8 

6.6 

87.4 

WRT 

Nominal 

0.7 

180 

0.0 

0.00 

0.00 

0.00 

0.00 

0.00 

0.0 

0.0 

0.0 

Nominal 

34 



3.0 

0.89 

- 0.85 

- 0.17 

- 0.20 

- 1.20 

- 1.3 

- 10.8 

46.9 



Nominal 

Runs 27&28 



Std Dev 

0.03 

0.02 

0.11 

0.13 

0.15 

0.2 

0.3 

3.7 

Condition 

Runs 27&28 



Average 

-7.69 

3.05 

1.55 

1.69 

5.61 

225.3 

172.3 

1011.7 


THICK values from Table 8 were compared with the ice tracings in Fig. lib, and this comparison revealed that 
quantitative changes agreed with the qualitative observations about ice feature change due to temperature variation. 
As expected, there was less variation in parameters at this freezing fraction of 0.7, versus what was observed at r| = 
0.3 in Table 7. Though the mass change was larger than in Table 7, it may still be within experimental accuracy 
associated with our mass measurement which had approximately 4 percent error associated with it. This accuracy 
estimate is based on cumulative experience from previous experiments where the authors have acquired ice shape 
mass measurements. Thus the value of 46.9 grams for a mass change would constitute about 4 percent of the ice 
mass acquired at the nominal condition (1011 gm). The value of 87.4 grams would constitute about an 8 percent 
change in ice mass. 


IV. Summary of Results 

Qualitative and quantitative results have been presented concerning the variation of LWC, MVD, and 
temperature and their effect on ice shape features. For the level of parameter variation selected for this study there 
were clearly discernible changes in ice shape features. But what if we had chosen a smaller level of variation in 
LWC, MVD, or temperature? How small could that variation be, before we would approach the level where ice 
shape feature changes would be on the order of the ice shape repeatability? This section will summarize the results 
from the previous sections in a way that attempts to answer these questions. 

Selected THICK data from the previous sections were subjected to an additional processing step. The icing limit 
data were excluded from this process because they were deemed to have too much variability due to possible errors 
inherent in the tracing process. Since errors in icing limits would also translate into errors in the calculated iced 
area, the iced area was excluded from this processing step. 

The process involved first normalizing the THICK calculated ice shape feature change. This was accomplished 
by dividing the calculated change value by the standard deviation for that ice feature. Next the parameter variation 
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was then divided by this normalized change value computed in the previous step. The end result was an indication 
of the parameter variation required to produce a standard deviation of change in the ice feature of interest. 

As an example: if the change in the upper icing horn were 16 degrees, and the standard deviation for this feature 
was 1.6, then the normalized change would be 10. If this change was caused by a -52pm variation in MVD (relative 
to the nominal condition), then the value of variation (-52) would be divided by the normalized change (10), to yield 
a result of -5.20. The interpretation of this result is that a decrease of 5.2 microns in drop size would produce the 
level of change in the ice feature (corresponding to a standard deviation). This processing step assumes that the 
change in ice feature is linear between the nominal condition and the variant condition. While the icing process is 
not linear, this may be acceptable as a first order approximation. The results from this process are summarized in 
Table 9 for variations in LWC, MVD, and temperature. 

Table 9 Parameter Variation to Produce Standard Deviation of Change in THICK Calculation 


Parameter 
Being Varied 

Description 

Run 

Parameter 

Change 

Lower 

Horn 

Thickness 

(inches) 

Upper 

Horn 

Thickness 

(inches) 

Upper 

Lower Horn Horn Angle 
Angle (deg) (deg) 

Mass (gm) 


Std Dev 



0.12 

0.06 

1.6 

1.6 

14.6 

LWC 

Change in 

20 

-0.10 

-0.14 

-0.23 

-7.7 

6.6 

-124.3 

(g/m 3 ) 

value 

21 

0.10 

0.01 

0.14 

11.5 

-5.4 

242.7 

n = o.3 

g/rrL per Std 

20 


0.084 

0.025 

0.020 

-0.024 

0.012 

V = 180 knots 

Dev 

21 


1.268 

0.042 

0.014 

-0.029 

0.006 


Std Dev 



0.11 

0.13 

0.2 

0.3 

3.7 

LWC 

Change in 

29 

-0.10 

-0.17 

-0.29 

2.7 

6.6 

-192.8 

(g/m 3 ) 

value 

30 

0.10 

0.12 

0.18 

-5.1 

-1.3 

412.6 

n = 0.7 

g/m 3 per Std 

29 


0.062 

0.044 

-0.007 

-0.004 

0.002 

V = 180 knots 

Dev 

30 


0.088 

0.071 

-0.004 

-0.022 

0.001 


Std Dev 



0.12 

0.06 

1.6 

1.6 

14.6 

MVD 

Change in 

22 

-52 

-0.13 

-0.02 

-15.4 

16.0 

-162.7 

(pm) 

value 

23 

50 

-0.22 

0.24 

11.8 

-3.8 

70.9 

n = o.3 

pm per Std 

22 


46.7 

132.7 

5.3 

-5.2 

4.7 

V = 180 knots 

Dev 

23 


-26.7 

12.2 

6.6 

-21.1 

10.3 


Std Dev 



0.11 

0.13 

0.2 

0.3 

3.7 

MVD 

Change in 

31 

-52 

-0.39 

-0.43 

-10.2 

14.3 

-327.6 

(pm) 

value 

32 

50 

0.27 

0.26 

-3.9 

-2.6 

263.7 

n = 0.7 

pm per Std 

31 


13.9 

15.6 

1.0 

-1.0 

0.6 

V = 180 knots 

Dev 

32 


19.5 

24.5 

-2.5 

-5.5 

0.7 


Std Dev 



0.12 

0.06 

1.6 

1.6 

14.6 

Ts 

Change in 

24 

-3.0 

0.10 

0.21 

-22.3 

24.9 

11.3 

(°F) 

value 

25 

3.0 

-0.65 

-0.17 

31.4 

-28.6 

-18.8 

CO 

o 

II 

c - 

°F per Std 

24 


-3.4 

-0.8 

0.2 

-0.2 

-3.9 

V = 180 knots 

Dev 

25 


-0.5 

-1.1 

0.1 

-0.2 

-2.3 


Std Dev 



0.11 

0.13 

0.2 

0.3 

3.7 

Ts 

Change in 

33 

-3.0 

0.10 

0.03 

-2.8 

6.6 

87.4 

(°F) 

value 

34 

3.0 

-0.17 

-0.20 

-1.3 

-10.8 

46.9 

n = o.7 

°F per Std 

33 


-3.1 

-13.1 

0.2 

-0.1 

-0.1 

V = 180 knots 

Dev 

34 


-1.9 

-2.0 

-0.4 

-0.1 

0.2 


The rows in Table 9 labeled "g/m 3 per Std Dev" contain the results of the processed THICK ice feature 
calculations for LWC variation. The values in these rows represent the amount of change in LWC required to 
produce a change in the THICK ice shape feature equal to the measured standard deviation of that feature. Values 
significantly less than the 0.1 g/m 3 LWC variation of this study would be very desirable, and could imply that 
smaller LWC content variations would yield discernable quantitative changes in these ice features. This same 
rationale applies to the MVD, and temperature variation data in this table. 
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Some of the calculated values have been highlighted if they had an absolute value less than 25 percent of the 
parameter variations selected for this study (0.1 g/m 3 , 50 um, and 3.0 °F). This was intended to identify potential 
ice features which would be sensitive enough to smaller parameter changes. The ice horn angles and mass seem to 
be two ice shape features which would be sensitive to smaller LWC, MVD, or temperature variations in this study. 

Table 9 represents a preliminary attempt to relate the sensitivity of quantitative ice shape features to variations in 
icing parameters like LWC, MVD, or temperature. There are recognized limitations in what can be ultimately 
concluded from the data presented in this report, but the results suggest quantitative measurements of ice shape 
features have potential for aiding in the definition of requirements, particularly for SLD ice protection codes. 

V. Conclusions 

A sensitivity study was conducted to evaluate the effect of LWC, MVD, and temperature changes on ice shape 
features. This study employed the following parameter variations: LWC = +/- 0.1 g/m 3 , MVD = +/- 50pm, and 
temperature = +/- 3 °F. Ice shape tracings, and ice shape mass measurements were acquired for each test condition. 
These data were then qualitatively and quantitatively analyzed. 

Preliminary results suggest that smaller parameter variations might produce discernible ice feature changes. 
These variations were found to be less than 0.1 g/m 3 for LWC, on the order of 10pm for MVD, and less than 1.0 °F 
for temperature. 

Additional conclusions drawn from this investigation are listed below. 

(1) The level of LWC, MVD, and temperature variation used in this study produced changes in ice shape 
features which were readily discernible using both qualitative and quantitative methods. The data suggest that some 
ice features like ice horn angle and mass may be sensitive enough to allow smaller LWC, MVD, and temperature 
variations to be used in future studies, but the study was not conclusive enough to establish the minimum parameter 
variation that could practically be achieved. The following uncertainties are believed to be factors that will 
practically limit the smallest achievable parameter variation: 

• Uncertainties due to small perturbations in parameters being held constant, but causing an undesired 
change in the ice shape of a variant parameter (e.g., MVD varied but slight changes occur in LWC, and 
temperature which are supposedly held constant). 

• Uncertainties in icing tunnel spray calibration definition of MVD and LWC values. 

• Uncertainties in LWC due to cloud uniformity effects. 

(2) Quantitative methods were successfully applied to the evaluation of ice shape features in this study, and were 
considered to be indispensable to the analysis process. It is believed that the quantitative methods are required to 
discern subtle changes in ice shape features, and facilitate comparisons over a range of different conditions. The use 
of quantitative methods removes the subjectivity that is typically associated with qualitative evaluations of ice shape 
feature changes. The THICK analysis program proved to be a very valuable analysis tool. 

(3) A simple method was used to estimate the level of LWC, MVD, and temperature variation associated with 
the repeatability error in selected THICK calculations (horn thickness, horn angles, and mass). It was found that 
changes in upper and lower horn angles, and mass seemed to have the greatest sensitivity to MVD, LWC, or 
temperature variations. 

(4) Accounting for the fact that LEWICE does not predict feather growth, there was generally good agreement 
between the trends in ice shape feature changes predicted by LEWICE, and what was observed from reviewing 
actual ice shape tracings. It was desired to include a quantitative evaluation of LEWICE predicted ice shapes for the 
conditions of this study, but that was not possible due to space limitations. 

(5) Errors in the ice shape generation and measurement process really need to be carefully controlled. Some 
suggested improvements to what was employed in this study are listed below: 

• Develop a better estimate of the repeatability of the nominal (or reference) conditions using more than 
two repeat runs. 

• The repeatability for the nominal condition was determined, but should probably have been done for 
variant conditions as well. 

• Some guidelines should be established concerning where to stop tracing, as you approach the icing 
limits 

• Develop some methodology to ensure the ice tracing template is properly positioned in the ice tracing 
cut on the airfoil 
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(6) It should be mentioned that the results of this study have limitations associated with only evaluating one 
airfoil. It is recognized that different airfoil characteristics (geometry, chord, etc) could exhibit differences relative 
to what has been reported here. 

Based on what has been learned from this study, there are a number of recommendations that can be made 
regarding further work in this area. 

• Complete the analysis of LEWICE predicted ice shapes with the THICK program, and compare the 
predicted ice feature changes with actual changes from this experimental study. 

• Assess the potential use of design of experiments principles for use in defining future testing. As part 
of this assessment, it may also be prudent to tiy and estimate the overall experimental uncertainty 
associated with ice feature changes (accounting for uncertainties in tunnel spray calibration, 
uniformity, ice shape tracing errors, ...etc.) 

• Develop some preliminary estimates of the needed resolution in ice feature change measurements, 
based on iced aerodynamic performance factors. It is envisioned this might be accomplished by 
reviewing the large body of work that exists on iced aerodynamic performance, and possibly 
complimented by conducting appropriate CFD investigations. These estimates could then be used to 
help define the desired level of parameter variations for future sensitivity tests. 

• Consider conducting additional sensitivity experiments on other airfoil geometries. 
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